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Abstract. In January 2013, North China Plain experienced 
several serious haze events. Cimel sunphotometer measure- 
ments at seven sites over rural, suburban and urban regions 
of North China Plain from 1 to 30 January 2013 were used 
to further our understanding of spatial-temporal variation 
of aerosol optical parameters and aerosol radiative forcing 
(ARF). It was found that Aerosol Optical Depth at 500 nm 
(AOD5oonm) during non-pollution periods at all stations was 
lower than 0.30 and increased significantly to greater than 
1.00 as pollution events developed. The Angstrom exponent 
(Alpha) was larger than 0.80 for all stations most of the time. 
AOD5oonm averages increased from north to south during 
both polluted and non-polluted periods on the three urban 
sites in Beijing. The fine mode AOD during pollution peri- 
ods is about a factor of 2.5 times larger than that during the 
non-pollution period at urban sites but a factor of 5.0 at sub- 
urban and rural sites. The fine mode fraction of AOD675nm 
was higher than 80% for all sites during January 2013. The 


absorption AOD 675 nm at rural sites was only about 0.01 dur- 
ing pollution periods, while ~ 0.03-0.07 and 0.01-0.03 dur- 
ing pollution and non-pollution periods at other sites, respec- 
tively. Single scattering albedo varied between 0.87 and 0.95 
during January 2013 over North China Plain. The size dis- 
tribution showed an obvious tri-peak pattern during the most 
serious period. The fine mode effective radius in the pollu- 
tion period was about 0.01-0.08 pm larger than during non- 
pollution periods, while the coarse mode radius in pollution 
periods was about 0.06-0.38 pm less than that during non- 
pollution periods. The total, fine and coarse mode particle 
volumes varied by about 0.06-0.34 pm 3 , 0.03-0.23 pm 3 , and 
0.03-0.10 pm 3 , respectively, throughout January 2013. Dur- 
ing the most intense period (1-16 January), ARF at the sur- 
face exceeded — 50 Wm -2 , — 180 Wm -2 , and —200 W m 2 
at rural, suburban, and urban sites, respectively. The ARF 
readings at the top of the atmosphere were approximately 
—30 W m 2 in rural and — 40-60 Wm -2 in urban areas. 
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Positive ARF at the top of the atmosphere at the Huimin sub- 
urban site was found to be different from others as a result of 
the high surface albedo due to snow cover. 


1 Introduction 

Aerosol particles can result in direct radiative forcing and 
they also have indirect effects on clouds (Hansen et al., 
1997). Despite a broad body of research on aerosols, their 
concentrations and optical properties remain one of the 
largest sources of uncertainty in current assessments and pre- 
dictions of global climate change (Hansen et al., 2000; Ra- 
manathan et al., 2001). Aerosol particles not only affect the 
climate, they can also cause other environmental and public 
health problems (Watson, 2005). 

China has experienced rapid economic growth over the 
past three decades, characterised by expanding industrializa- 
tion, explosive growth in cars and frequent anthropogenic ac- 
tivities. These changes have resulted in high aerosol loading 
levels (Qiu and Yang, 2000; Luo et al., 2001; Zhang et al., 
2013b; Li et al., 2013), decreases in visibility (Che et al., 
2007), solar radiation (Che et al., 2005; Liang and Xia, 2006; 
Xia, 2010) and air quality (Xu et al., 2002), especially in ar- 
eas of eastern China such as the Jing-Jin-Ji (Beijing, Tianjin 
and Hebei) region (Xu et al., 2005; Li et al., 2007a, b, 2011; 
Ding et al., 2008), the Yangtze River Delta and parts of south- 
ern China such as around the Pearl River Delta (Wang et al., 
2003; Gao et al., 2011). 

Employing ground-based measurement networks is very 
useful and accurate when investigating aerosol optical prop- 
erties (Holben et al., 2001). There are several well-known 
ground-based aerosol-monitoring networks that use sunpho- 
tometers such as the Aerosol Robotic Network- AERONET 
(Holben et al., 1998), PHOtometrie pour le Traitement 
Operationnel de Normalisation Satellitaire (PHOTONS) 
(Goloub et al., 2008), AEROCAN (Bokoye et al., 2001), 
China Aerosol Remote Sensing Network (CARSNET) (Che 
et al., 2009a), and SKYNET network (Takamura and Naka- 
jima, 2004). These networks have been used to measure di- 
rect and diffuse solar radiation and derived aerosol optical 
properties for the purpose of aerosol radiative forcing stud- 
ies (Dubovik et al., 2002; Kim et al., 2004). There has been a 
wealth of research on the optical properties of aerosols above 
the Jing-Jin-Ji region of China in recent years (Cheng et al., 
2008; Che et al., 2008; Fan et al., 2009; Chen et al., 2013; Li 
et al., 2013; Xia et al., 2013), which has resulted in a com- 
prehensive understanding of the aerosol optical properties of 
this region (Eck et al., 2005). These researches will help not 
only to uncover the regional climate effect of aerosol parti- 
cles and the mechanism of pollution event formation, but also 
will allow managers to control air quality better in the future. 

The most frequent and serious particle matter (PM) 
pollution of this century appeared in January 2013 
in eastern China, especially in the densely populated 


North China Plain, which caused considerable pub- 
lic concerns (http://edition.cnn.com/2013/01/19/world/asia/ 
china- florcruz- pollution). The greatest instantaneous con- 
centration of PM2.5 reached 1000 pgm -3 in some heavily 
polluted areas of Beijing (Zhang et al., 2013a). The objec- 
tive of this study is to investigate aerosol optical properties 
and their direct radiative forcing in rural, suburban and urban 
regions of North China Plain during this serious pollution 
episode. These high temporal resolution aerosol optical prop- 
erty data could be used in the validation of haze-fog mod- 
elling in future. The main findings include: (1) temporal vari- 
ations in aerosol optical properties over North China Plain; 
(2) inter-comparison of aerosol optical properties between 
pollution and non-pollution periods in North China Plain; 
and (3) the direct radiative forcing characteristics caused by 
pollution aerosol particles in North China Plain. 


2 Site description, instrumentation and 
analytical method 

2.1 Site description and instrumentation 

Seven sites in North China Plain providing sun and sky scan- 
ning radiometer (Cimel Electronique CE-318) observations 
were selected in this study (Fig. 1). The geographical infor- 
mation and site description are shown in Table 1 . The Shang- 
dianzi rural site is located 150 km to the northeast of Beijing, 
which is regarded as the regional atmospheric background 
station because it is far away from large anthropogenic emis- 
sion sources. Therefore, the measurements from Shangdi- 
anzi represent the basic characteristics of aerosols in North 
China Plain, especially in the Beijing-Tianj in-Hebei (Jing- 
Jin-Ji) Metropolitan region (Tang et al., 2007; Che et al., 
2009b). Stations such as the Institute of Remote Sensing 
and Digital Earth (RADI), the Chinese Academy of Mete- 
orological Sciences (CAMS), and the Beijing Meteorologi- 
cal Observatory (Nanjiao) are located in the Beijing urban 
area. Measurements at these three stations not only reflect 
the aerosol characteristics in urban North China Plain, but 
also represent the aerosol characteristics in the north, mid- 
dle and southern parts of the greater Beijing urban area, re- 
spectively (Che et al., 2009a). Xianghe and Huimin are two 
suburban sites in North China Plain, the measurements taken 
there reflecting the aerosol characteristics of suburban North 
China Plain (Xia et al., 2005). Similar to RADI, CAMS, 
and Nanjiao stations, the Tianjin site is located in the cen- 
tre of Tianjin city, so that these measurements also reflect the 
aerosol characteristics of urban North China Plain. Among 
these seven sites, RADI, CAMS, and Xianghe stations are the 
AErosol RObotic NETwork- AERONET stations (Holben et 
al., 1998; Li et al., 2007, 2013). Shangdianzi, CAMS, Nan- 
jiao, Tianjin, and Huimin belong to the China Aerosol Re- 
search Network-CARSNET (Che et al., 2009a). The CAMS 
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Fig. 1. Distribution of the seven sun photometer sites in North China 
Plain. Wind rose maps were plotted at Shangdianzi, Nanjiao, Tian- 
jin, and Huimin, using observation data on the hourly meteorologi- 
cal wind speed and direction. 


is the main station of CARSNET but the measurements are 
also uploaded to the AERONET data archive. 

The Cimel Electronique CE-318 sunphotometer used in 
this research has eight channels, including 1640, 1020, 870, 
670, 500, 440, 380 and 340 nm and a 940 nm water vapour 
channel with a 1.2° full field-of-view (Holben et al., 1998). 
Measurements at 1640, 1020, 870, 670, 500, 440, 380 and 
340 nm can be used to retrieve AOD, and measurements at 
940 nm quantify in centimetres the total precipitable water 
vapour. The raw AOD is cloud-screened according to the 
method of Smirnov et al. (2000). The total uncertainty asso- 
ciated with the aerosol optical depth is approximately 0 . 01 - 
0.02 (Eck et al., 1999). The sun photometer can also de- 
tect sky radiance for solar almucantar scenario or principal 
plane scenario measurements within a 1 . 2 ° full field-of-view 
at 440, 675, 870, and 1020 nm. These sky radiance data are 
measured at 24 pre-defined scattering angles at regular time 
intervals, which are used to retrieve aerosol optical param- 
eters such as single scattering albedo, size distribution, and 
refractive index, etc. (Dubovik and King, 2000; Dubovik et 
al., 2006). 

In order to retain the accuracy of aerosol optical re- 
trieval parameters, AERONET undertakes the Langley cali- 
bration at Mauna Loa (MLO, 19.54° N, 55.58° W, 3397.0 m) 
and inter-comparison calibration and sphere calibration is 
conducted at the Goddard Space Flight Center (GSFC, 
38.99° N, 76.84° W, 87.0 m) periodically (Holben et al., 
1998). AERONET-PHOTONS compiles the Langley calibra- 
tion at the Izana Observatory in the Canary Islands (INM, 
Spain, 28.31° N, 16.50° W, 2391.0m) and Lille University 
(France, 50.61° N, 3.14° E, 60.0 m) is responsible for inter- 
comparison calibration and sphere calibration (Gloub et al., 
2007). CARSNET records the Langley calibration at the 


Izana Observatory (INM, Spain), following the calibration 
protocol used by NASA. The inter-comparison calibration 
is carried out at CAMS using an inter-comparison method 
(Che et al., 2009a). The sphere calibration is performed by 
the CARSENT sphere every year to ensure the accuracy of 
the sky irradiance measurement (Tao et al., 2014). 

2.2 Data description and analytical method 

In this study, the meteorological data from 1 to 30 January 
2013, including hourly relative humidity, wind speed and 
direction, has been obtained from the China Meteorologi- 
cal Administration (CM A). The hourly PM 2.5 data for the 
same period in Beijing were supplied by the China National 
Environmental Monitoring Centre (http://www.cnemc.cn/). 
There are about 30 stations located in Beijing. The PM 2.5 
mass concentration data in this study was the average of 
all these stations, which could be representative of the en- 
tire Beijing. The wind frequency map was calculated us- 
ing the 10 min wind data at Shangdianzi, Nanjiao, Tianjin 
and Huimin (Fig. 1). The daily averaged PM 2 . 5 , relative hu- 
midity and wind speed and wind direction at Beijing were 
calculated and are shown in Fig. 2. The Cimel sun pho- 
tometer measurement was also selected for the same Jan- 
uary period. Aerosol optical property and its radiative forcing 
data of RADI, CAMS, and Xianghe are from Level 1.5 data 
(cloud- screened) in the AERONET website (http://aeronet. 
gsfc.nasa.gov/). Aerosol optical property and its radiative 
forcing data for Shangdianzi, Nanjiao, Tianjin and Huimin 
collected from CARSNET were processed in a method sim- 
ilar to that used with AERONET (see next section). 

2.3 Retrieval methods of Cimel measurement 
of CARSNET 

The AOD data of the four CARSNET sites at Shangdianzi, 
Nanjiao, Tianjin and Huimin were processed using the AST- 
Pwin software offered by Cimel Ltd. Co (Che et al., 2009a). 
This software is used to calculate Level 1.0 AOD (raw result 
without cloud- screening), Level 1.5 AOD (cloud- screened 
AOD based on the work of Smirnov et al., 2000) and the 
Angstrom Exponent (Alpha) between 440 nm and 870 nm. 
Che et al. (2009a) pointed out that there is good consis- 
tency between the CARSNET calculation results and those 
of AERONET/PHOTONS. The daily averaged AOD is less 
than 2.05 %, 0.46 %, 1.50 %, and 2.61 % at 1020 nm, 870 nm, 
670 nm, and 440 nm, respectively, based on comparison of 
the measurements between the three CARSNET instruments 
and AERONET/PHOTONS. 

The retrieval processing method of sky radiance measure- 
ments used by CARSNET in this study is very similar to the 
Version 2 AERONET algorithm (Dubovik et al., 2006). In 
the AERONET Version 2 algorithm, the Cimel sky radiance 
measurements in the almucantar plane (fixed elevation angle 
equal to solar elevation and a full 360° azimuthal sweep) at 
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Table 1 . Site location and description. 


Station name 

Lat (°) 

Long (°) 

Site description 

Shangdianzi 

40.65° N 

117.12° E 

Rural station, 393.0 ma.s.l., 150 km to the northeast of Beijing city 

RADI 

40.00° N 

116.38° E 

Urban station, 59.0 ma.s.l., in the north part of Beijing 

CAMS 

39.93° N 

116.32° E 

Urban station, 105.0 m a.s.l., in the central part of Beijing 

Nanjiao 

39.80° N 

116.47° E 

Urban station, 31.3 ma.s.l., in the south part of Beijing 

Xianghe 

39.75° N 

116.96° E 

Suburban station, 36.0 m a.s.l., 50 km to the east of Beijing 

Tianjin 

39.10° N 

117.17° E 

Urban station, 3.3 ma.s.l., in the center of Tianjin City, 120 km to the southeast of Beijing 

Huimin 

37.48° N 

117.53° E 

Suburban station, 1 1 .7 m a.s.l., 280 km to the southeast of Beijing 


440, 675, 870, and 1020 nm (nominal wavelengths) in con- 
junction with the direct sun measured AOD at these same 
wavelengths (AOD 44 o n m should be larger than 0.40) were 
used to retrieve column-integrated aerosol size distributions 
(dV(r)/dln(r) from 0.05 to 15 pm) and refractive indices, sin- 
gle scattering albedo, retrieved AOD, effective and median 
radii, volume concentration of fine and coarse mode aerosols 
at low solar elevations (solar zenith angles between 50 and 
80 degrees) (Eck et al., 2010). 

The radiative forcing and the forcing efficiency data at 
all sites in this study were calculated by the same radiative 
transfer module which has been integrated into the opera- 
tional AERONET inversion code (Garcia et al., 2012). This 
module uses the detailed size distribution, complex refrac- 
tive index and fraction of spherical particles retrieved from 
Cimel sky radiance measurements in the almucantar plane 
(Dubovik and King, 2000; Dubovik et al., 2006) and pro- 
vides the fluxes and aerosol radiative forcing (ARF) values 
using the Discrete Ordinates DISORT approach (Nakajima 
and Tanaka, 1988; Stamnes et al., 1988). Note that instan- 
taneous ARF for solar zenith angle of 50-80° is presented 
here. 


3 Results and discussion 

3.1 Analysis of PM 2.5 and meteorological data 

Before analysing the aerosol optical properties and its 
radiative forcing, we take measurements at Beijing as 
an example to investigate the temporal variation of 
PM 2.5 and meteorological elements because these two 
factors contribute substantially to air pollution. Figure 2 
shows the daily mean variation of PM 2 . 5 , relative hu- 
midity, and wind speed and direction in Beijing. Ac- 
cording to the ambient air quality standards of China 
(GB 3095-2012: http://kjs.mep.gov.cn/hjbhbz/bzwb/dqhjbh/ 
dqhj zlbz/20 1 203/t20 1 20302_224 1 65 .htm), which will be 
implemented from 1 January 2016, air quality is regarded as 
a “pollution day” if the 24 h average of PM 2.5 mass concen- 
tration exceeds 75pgm -3 . Thus, applying these criteria to 
the month of January 2013, there are 21 pollution days (pol- 
luted period) and only 10 non-pollution days (non-polluted 


period) over Beijing (Fig. 2a). The most intense pollution 
episode occurred from January 11th to 17th with the maxi- 
mum PM 2.5 mass concentration > 500 pg m -3 on 13 January. 
Temporal variation of relative humidity was consistent with 
that of PM 2.5 (correlation coefficient being 0.58). On the 
pollution days, the relative humidity usually exceeded 50 % 
(Fig. 2b). Hennigan et al. (2008) pointed out that the for- 

o 

mation of secondary aerosol species such as N0 3 , S0 4 
and secondary organic compounds could be enhanced under 
high-humidity conditions similar to the in-cloud process. In 
general, temporal variation of wind speed showed a negative 
relation with the PM 2.5 (correlation coefficient being —0.25). 
During the polluted period, the wind speed was usually less 
than 2.0ms _1 (Fig. 2c). However, a higher negative relation 
(correlation coefficient being —0.43) between wind speed 
and the PM 2.5 on the next day also occurred in January 2013 
in Beijing, indicating that wind speed is an important factor 
affecting PM 2 . 5 . Both low and high wind speed affected the 
accumulation and diffusion of aerosol particles during Jan- 
uary 2013. The temporal wind direction variation is shown 
in Fig. 2d. Fow PM 2.5 concentrations usually coincide with 
northerly winds (0-90°, 270-360°) on the previous day (e.g. 
1-4, 6 , 9 and 20 January). The wind from the clean area to the 
north can reduce the pollutants and transport them to down- 
stream region while, during the pollution period, wind direc- 
tion on the previous day was close to southerly where an- 
thropogenic activity is concentrated and frequent, especially 
in Hebei Province (Xia et al., 2007a). Apart from local emis- 
sions in Beijing, additional aerosol particles transported to 
Beijing can cause deterioration of the air quality. On the wind 
rose map (Fig. 1), the majority of wind frequencies during 
January 2013 at Beijing, Tianjin and Huimin come from the 
inland North China Plain where there are many industrial ac- 
tivities. Shangdianzi station is different from the others be- 
cause the majority of wind frequencies are from the eastern 
mountain regions and the southwestern urban regions. 

3.2 Spatial variation of AOD, Angstrom exponent 
(Alpha), and volume size distribution (dV/dlnr) 
during January 2013 in North China Plain 

As the spatial and temporal variation of AOD at 500 nm can 
be seen on Fig. 3, AOD variation patterns were consistent 
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Date (January 2013) 


Fig. 2. Daily averaged PM 2.5 (the dashed line indicates 75 pg m -3 ) 
(a), relative humidity (b), wind speed (c), wind direction (the upper 
red dash line indicates the wind from west and the bottom one indi- 
cates the wind from east) (d) at Beijing. The shaded area indicates 
daily PM2.5 exceeding 75 pgm -3 . 


in all rural, suburban, and urban regions during the pollution 
month of 2013. AOD values in the non-pollution period at 
all stations were lower than 0.30 and increased significantly 
to more than 1.00 as the pollution events occurred. This con- 
sistent variation indicates that intense pollution is not a local 
but a regional phenomenon in North China Plain. During the 
whole of January 2013, the daily AODsoonm varied by about 
0.15-2.60 (minimum 0.09, maximum 3.60) at CAMS, 0.13- 
2.79 (minimum 0.09, maximum 3.25) at Nanjiao, 0.12-2.70 


(minimum 0.08, maximum 1.86) at RADI, 0.06-2.00 (mini- 
mum 0.05, maximum 2.51) at Shangdianzi, 0.08-2.12 (mini- 
mum 0.06, maximum 2.12) at Xianghe, 0.17-3.17 (minimum 
0.06, maximum 2.12) at Tianjin, and 0.17-2.39 (minimum 
0.10, maximum 2.87) at Huimin, respectively. AOD varia- 
tion at the Shangdianzi station was consistent with those at 
RADI, CAMS, and Nanjiao, but the values were systemat- 
ically lower than that in these three urban sites, suggesting 
that there are fewer anthropogenic pollution sources in this 
region. For the three urban sites of RADI, CAMS and Nan- 
jiao located in Beijing, the AOD during the intense periods 
from 10 to 16 January and 26 to 29 January, AOD reached 
a level as high as 2.00-3.00, which is similar to the previous 
results. Li et al. (2013) pointed out that the averaged AOD 
at 440 nm could reach 3.2 during thick haze events in Bei- 
jing. AOD variations at all four of the Beijing sites (Shangdi- 
anzi, RADI, CAMS, and Nanjiao) are consistent with those 
of PM 2.5 (Fig. 2a). Li et al. (2013) pointed out that the corre- 
lation between AODsoonm and PM 2.5 could reach 0.93 dur- 
ing haze events. In this study, the correlation between daily 
AOD 5 oonm and PM 2.5 during January 2013 was 0.82, 0.75 
and 0.70 for RADI, CAMS, and Nanjiao, respectively, sug- 
gesting that the low level of atmosphere is well mixed during 
the pollution episode. 

Figure 4 shows the spatial variation of Alpha during Jan- 
uary 2013. It is obvious that the Alpha was larger than 0.80 
for all stations and for most of the time, which suggests that 
the small aerosol particles was dominant during the pollution 
month over North China Plain. This result is similar to some 
previous ones (Eck et al., 2005; Xia, 2007a). During the two 
periods 12-16 and 24-29 January, the Alpha increased grad- 
ually at most stations, which suggests that the aerosol size 
became smaller as the pollution process developed. On 24 
January, the Alpha values at Shangdianzi, RADI, CAMS, and 
Nanjiao were all lower than 0.80, which suggests that coarse 
particles were dominant. This could be due to the decreasing 
relative humidity and increasing wind speed (Fig. 1). Coarse 
particles could be from more natural and anthropogenic dust 
emission with increasing wind speed. The PM 2.5 decreased 
from ~ 320 pg m -3 to ~ 120 pg m -3 on 23-24 January. 

Figure 5 expresses the size distribution and spatial varia- 
tion of particles at the three urban sites in Beijing. On pol- 
luted days, fine mode particles were clearly higher than in 
the non-polluted ones. In general, both fine and coarse mode 
particle volume increased gradually from 2 to 6 January at 
the three Beijing urban sites. In contrast to other polluted pe- 
riods with only two peaks at the radius about 0 . 10 - 0.20 pm 
and 1.20-1.40 pm, the size distribution showed an obvious 
tri-modal pattern during the most intense period (10-16 Jan- 
uary) at the Beijing urban area. There were three peaks at 
the radius about 0.15 pm, 0.30-0.50 pm, and 1.10-1.20 pm. 
This interesting phenomenon probably reflects the hygro- 
scopic characteristics of fine particles. The relative humid- 
ity was more than 60% from 10 to 16 January. According 
to the results by Zhang et al. (2013a), the mean measured 
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Date (January 201 3) 


Fig. 3. Daily averaged AOD (green) and the 5th and 95th percentile 
box plots for the month of January 2013 at seven sites in North 
China Plain. 


non-refractory submicron particle mass concentration at the 
Institute of Atmospheric Physics in Beijing (midway be- 
tween the RADI and CAMS in this study) was composed 
of organics (49.8 %), sulphate (21.4 %), nitrate (14.6 %), am- 
monium (10.4%), and chloride (3.8%). These hygroscopic 
compositions could be one of the major reasons for the tri- 
modal size distribution during the most intense haze period 
from 10 to 16 January. Another possible reason of the bi- 
mode submicron size distribution could be related to the 
haze-fog and cloud processing (Li et al., 2011; Eck et al., 
2012). The smaller mode (0.12 to 0.20 micron meter) may 
be interstitial aerosol that were not modified by incorpora- 
tion in droplets and/or aerosol that are less hygroscopic in 
nature. The larger one is due to the haze-fog dissipation. 



Date (January 201 3) 

Fig. 4. Daily averaged Alpha (green) and its 5th and 95th percentile 
box plots for January 2013 at seven sites in North China Plain. 


3.3 Inter-comparison of aerosol optical properties 
during pollution and non-pollution periods 

According to the PM2.5 spatial variation in Beijing, the Jan- 
uary 2013 period is divisible into the two periods, i.e. non- 
pollution and pollution period. In this section, aerosol opti- 
cal characteristics between non-pollution and pollution pe- 
riods are compared. Figure 6 shows the inter-comparison 
of AOD5oonm (a), Alpha (b), fine mode AOD675nm (c), the 
fine mode AOD675 nm fraction (d), absorption AOD675 nm (e), 
and SSA675 nm (f) during the two periods. AOD5oonm dur- 
ing the pollution period is a factor of 2-4.3 times the non- 
pollution period. During the polluted period, the AODsoonm 
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Fig. 5. Waterfall plots of daily aerosol volume size distribution dur- 
ing January 2013 at the three urban sites of RADI, CAMS and Nan- 
jiao, located in Beijing. 


averaged 0.45 zb 0.48, 0.66 zb 0.55, 0.98 zb 0.81, 1.36 zb 0.84, 
1.15 zb 0.66, 1.43 zb 0.90, and 1.21 zb 0.70 at Shangdianzi, 
RADI, CAMS, Nanjiao, Xianghe, Tianjin, and Huimin, 
respectively, and 0.12 zb 0.10, 0.28 zb 0.16, 0.35 zb 0.18, 
0.39 zb 0.25, 0.26 zb 0.16, 0.42 zb 0.25, and 0.59 zb 0.48, re- 
spectively, during the non-polluted period. The two urban 
sites of Nanjiao (1.36 zb 0.84) and Tianjin (1.43 zb 0.90) dur- 
ing the polluted period was a little higher than at the subur- 
ban sites of Xianghe (1.15 zb 0.66) and Huimin (1.21 zb 0.70), 
which is probably due to more local emissions at the urban as 
compared to that at the suburban sites. As for the three urban 
sites of RADI, CAMS and Nanjiao in Beijing, AODsoonm av- 


erages increased from the north (RADI) to the southern part 
(Nanjiao) during both the pollution and the non-pollution pe- 
riods, which clearly reflected the effect of local emissions. 
The RADI site is located outside of Beijing’s 4th ring road 
and near to the Olympic Forest Park where few industrial 
emissions occur. The CAMS is located between the 2nd and 
3rd Beijing ring road with heavier traffic loads compared 
to the RADI site, while the Nanjiao site is at a height of 
31ma.s.l. and only 200 m from the 5th Beijing ring road, 
where there are many factories near to the major ring road. 
Thus the higher local pollution may well contribute to the 
higher AOD found in the northern sites of CAMS and RADI. 

The Alpha readings during both the pollution and non- 
pollution periods are greater than 1.00 at all 7 sites, which 
suggest that the small particles are dominant. Alpha at 
Shangdianzi (1.06 zb 0.15) is lower than at the other sites, 
which indicates that the aerosol size in the rural site is larger 
than in the urban and suburban sites. This is because Shang- 
dianzi is located about 100 km northeast of urban Beijing and 
surrounded by a mountain chain where is sparsely populated 
(Hanel et al., 2012). Thus fewer fine anthropogenic aerosol 
particles could affect this area. In contrast to the AOD, Al- 
pha decreased from the north (RADI) to the southern location 
(Nanjiao) during both the polluted and non-polluted periods 
in Beijing, suggesting that larger size aerosol particles exist 
in the southern suburbs of Beijing. 

The fine mode fraction of AOD was computed based on 
the retrieved size distributions and spectral refractive indices 
after the Dubovik and King (2000) algorithm applied to al- 
mucantar scans, assuming bimodal size distribution. This pa- 
rameter is defined as the fraction of the total optical depth 
attributed to the fine mode of the aerosol size distribution, 
which could be used to discriminate aerosols with similar so- 
lar extinction but different sizes. The maximum radius of the 
fine mode volume distribution was defined as the minimum 
between the fine and coarse modes in the retrieved size distri- 
bution, over the defined range limits of 0.44-0.99 m radius. 
The minimum radius for the fine mode is fixed at a 0.05 mm, 
which is the minimum radius of the almucantar inversion 
(Eck et al., 2010). During January 2013, large fine mode 
AODs were observed in Tianjin (0.86 zb 0.31) and Huimin 
(0.48 zb 0.36) (Fig. 6b). Xianghe was similar to Nanjiao with 
a fine mode AOD of 0.35 zb 0.29. The fine mode AOD at 
Shangdianzi during the same month was just 0.14. The fine 
mode AOD in the pollution period at the urban sites of RADI, 
CAMS and Nanjiao was about a factor of 2.5 times the count 
in the non-pollution period, and a factor of about 5 on the 
suburban and rural sites. It can be seen in Fig. 6d that the 
fine mode fraction of AOD 675 nm (FMF) in January 2013 was 
larger than 80 % for all sites, and was even greater (> 90 %) 
at Tianjin and Huimin, indicating that the fine mode particles 
are the dominant contributor to the atmospheric extinction 
during January 2013 over North China Plain. This is similar 
to the results of Eck et al. (2010) and Xia et al. (2013), who 
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Fig. 6. Means of aerosol optical properties during the total period 
(January 2013), pollution period (daily PM 2.5 > 75pgm -3 ), and 
non-pollution period (daily PM 2.5 < 75pgm“ 3 ) at each site, (a) 
AOD at 500 nm, (b) Alpha between 440-870 nm, (c) Fine mode 
AOD at 675 nm, (d) Fine mode AOD fraction at 675 nm (Fine 
mode AOD 675 nm /Total mode AOD 675 nm ), (e) Absorption AOD at 
675 nm, (f) Single scattering albedo at 675 nm (SSA 675 nm ). 


pointed out that FMF could reach about 0.93 during most of 
the predominantly fine mode-dominated pollution cases. 

Figure 6e shows the absorption AOD 675 nm (AAOD) at all 
sites during both the pollution and non-pollution periods. The 
AAOD at Shangdianzi was only about 0.01 during the pollu- 
tion period while, for other sites, AAOD was about 0.03- 
0.07 and 0.01-0.03 for pollution and non-pollution periods, 
respectively. This indicates that there were few absorption 
aerosol particles at the rural site of Shangdianzi and more 


absorption aerosol particles at both suburban and urban sites 
during the pollution period. These absorption aerosol parti- 
cles may consist of black carbon and brown carbon emitted 
by coal combustion and vehicle exhausts (Zhang et al., 2012; 
Li et al., 2010). Li et al. (2013) found that the proportion of 
black carbon and brown carbon can reach between 25 and 
38 % during pollution events. 

The single scattering albedo (SSA) varied between 0.87 
and 0.95 in January 2013 at all seven sites (Fig. 6f). For 
the three Beijing urban sites and the two suburban sites of 
Xianghe and Huimin, the SSA values for the pollution pe- 
riod (0.88-0.93) are higher than for the non-pollution period 
(0.85-0.90), which could be due to the presence of water- 
soluble aerosols at higher relative humidity, resulting in fine 
mode particle growth and an increase in the light scattering 
coefficient (Kotchenruther and Hobbs, 1998). 

Figure 7 shows the total, fine and coarse mode effec- 
tive radii and volumes of aerosol particles during the pol- 
lution and non-pollution periods at all seven sites. The to- 
tal mode effective radii varies from about 0.28 to 0.36 pm, 
0.30 to 0.40 pm, and 0.23 to 0.41 pm for the whole period, 
the pollution period and the non-pollution period, respec- 
tively (Fig. 7a). The fine mode effective radii varies in the 
ranges 0.14-0.23 pm, 0.15-0.23 pm, and 0.13-0.15 pm for 
the whole period, the pollution period and the non-pollution 
period, respectively (Fig. 7b), and the values of the coarse 
mode effective radii are about 2.21-2.46 pm, 2.09-2.39 pm, 
and 2.10-2.62 pm for the same series of periods (Fig. 7c). 
The fine mode effective radius of the pollution period is about 
0.01-0.08 pm larger than that found in the non-pollution 
period, which reflects the fine mode particle hygroscopic 
growth mentioned above. In contrast, the coarse mode ra- 
dius in the pollution period was about 0.06-0.38 pm less than 
that in the non-pollution period. It could be speculated that 
the higher wind speed during non-pollution periods results in 
more coarse particles in the atmosphere such as, for instance, 
the fugitive dust (Wang et al., 2009) and fly ash emissions 
from coal burning (Yang et al., 2009). 

The total, fine and coarse mode particle volumes varied by 
0.06-0.34 pm 3 , 0.03-0.23 pm 3 , and 0.03-0.10 pm 3 , respec- 
tively, during the whole month of January 2013 (Fig. 7d- 
e). During the pollution period, particle volumes were about 
0.07-0.17 pm 3 , 0.03-0.12 pm 3 , and 0.01-0.05 pm 3 larger 
than in the non-pollution period, which is suggestive of more 
aerosol loading (both fine and coarse mode particles) of the 
atmosphere during pollution periods over North China Plain. 

3.4 Analysis of aerosol radiative forcing (ARF) and its 
efficiency 

As pointed by Garcia et al. (2012), the aerosol radiative forc- 
ing (ARF) provides the actual or total radiative effect of at- 
mospheric aerosols. To make a consistent comparison among 
AFR of atmospheric aerosols the aerosol radiative forcing ef- 
ficiency is a magnitude more appropriate. This parameter is 
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Fig. 7. Means of aerosol size distribution properties during the total 
period (January 2013), pollution period (daily PM 2.5 > 75 pg m -3 ), 
and non-pollution period (daily PM 2.5 < 75pgm -3 ) at each site, 
(a) Retrieved effective radius of total mode; (b) retrieved effective 
radius of fine mode; (c) retrieved effective radius of coarse mode; 
(d) retrieved aerosol particle volume of total mode; (e) retrieved 
aerosol particle volume of fine mode; and (f) retrieved aerosol par- 
ticle volume of coarse mode. 


defined as the rate at which the atmosphere is forced per unit 
of AOD. 

Figures 8 and 9 illustrate the spatial variations of ARF 
at the surface (ARF-BOA) and the top of the atmosphere 
(ARF-TOA) at the Shangdianzi, RADI, CAMS, Nanjiao 
and Huimin stations. As pollution levels became more in- 
tense during the periods 9-13 and 25-30 January, the ARF- 
BOA values progressively increased (Fig. 8). During the 
most intense period (10-16 January), the ARF-BOA ex- 
ceeded — 50 Wm -2 at Shangdianzi; — 200 Wm -2 at the 
three Beijing urban sites of RADI, CAMS, and Nanjiao; and 
— 180Wm -2 at the Huimin suburban site. The ARF-TOA 


varied with the same trends as the ARF-BOA at all sites dur- 
ing January 2013. During the most intense polluting period, 
the ARF-TOA values were — 30 Wm -2 at Shangdianzi, and 
about 40-60 Wm -2 at the three Beijing urban sites. How- 
ever, the AFR-TOA spatial variation at Huimin was very dif- 
ferent from the other sites. Before 4 January, the Huimin 
AFR-TOA was negative; this suggests that the aerosol par- 
ticles imposed a cooling effect at the top of the atmosphere. 
From 5 January, the ARF-TOA at the same site became pos- 
itive, which may have been a result of a heating effect by 
the aerosol particles at the top of the atmosphere. The rea- 
son for positive AFR-TOA is the presence of a large sur- 
face albedo (SA). (In this study, the surface albedo refers 
to records obtained from the MODerate-resolution Imaging 
Spectroradiometer (MODIS) albedo product (MCD43C3).) 
The precipitation records indicate that, from 20 December 
2012, there were 3 days with a snow cover at Huimin. Thus, 
the MODIS albedo product showed surface albedo at 440 nm 
for Huimin, varying from 0.3498 to 0.5238 throughout the 
whole month of January 2013, which is much larger than the 
data from CAMS of 0.0645 to 0.094 on the same dates. Gar- 
cia et al. (2012) points out that the effect of surface albedo 
on the AFR is critical for SA >0.30. The high surface albedo 
could have caused more shortwave radiation reflecting back 
to the atmosphere and heating it by the absorption effect of 
aerosol particles (Wang et al., 2011). Finally, the high sur- 
face albedo was the cause of the positive ARF at the top of 
the atmosphere. 

Figure 10 shows the ARF-BOA, the ARF-TOA and the ef- 
ficiency of both during the whole period and the pollution 
and non-pollution periods. The ARF-BOA varied at all seven 
sites across North China Plain from —32 to — 144 Wm -2 , 
—43 to — 144 Wm -2 , and —9 to —63 W m -2 for the whole, 
pollution and non-pollution periods, respectively, which in- 
dicates aerosol particles produce a cooling effect at this re- 
gion. The ARF-TOA varied at all seven sites in North China 
Plain from 11 to — 54 Wm -2 , 18 to —54 W m” 2 , and —3 to 
— 21 Wm -2 for the whole, pollution and non-pollution pe- 
riods, respectively. The ARF-TOA values are smaller than 
those in the ARF-BOA which suggests the aerosol parti- 
cles higher absorption reduces the solar energy available to 
be backscattered to TOA and keep a larger percentage in 
the atmosphere. The ARF-BOA efficiencies are more than 
— 150 W m -2 for all sites except Tianjin (— 133 zb 17 W m -2 ) 
during the whole, pollution and non-pollution periods. 
The largest ARF-BOA efficiencies occurred in RADI sites 
with values of — 240 ± 51 Wm -2 , — 226±54Wm -2 and 
— 256±43Wm -2 during the whole, pollution and non- 
pollution periods, respectively. This could be due to the 
lower SSA values (less than 0.88) found at RADI. More 
absorbing aerosols can induce larger ARF efficiency at sur- 
face (Garcia et al., 2012). The ARF-TOA efficiencies RADI, 
CAMS, Nanjiao and Xianghe show a set of similar lev- 
els as follows: —64 to — 87 Wm -2 , —55 to — 77Wm -2 
and —67 to — lOOWm -2 during the whole, pollution and 
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Date (January 2013) 

Fig. 8. Aerosol radiative forcing at the surface (ARF-BOA) and its 
5th and 95th percentile box plots during January 2013 at five sites 
in North China Plain. 


non-pollution periods, respectively. The Shangdianzi ARF- 
TOA efficiencies are the highest among all 7 stations 
with values of -121±21Wm- 2 , -114±23Wm~ 2 , and 
— 134±10Wm -2 during the whole, pollution and non- 
pollution period, respectively. This is because the minimum 
AOD was observed there and ARF decreases with AOD (Xia 
et al., 2007b). Garcia et al. (2012) pointed it clearly that even 
not significant uncertainties in the absolute values of radia- 
tive forcing at the low AOD can present more significant 
impact on relative to AOD values of forcing efficiency. The 
ARF-TOA efficiency at Huimin during the pollution periods 
was 18.89 ± 18.85 Wm“ 2 , but -32.36 ± 40.78 Wm“ 2 dur- 
ing the non-pollution periods. This positive ARF-TOA effi- 
ciency under the polluted condition is mainly due to the large 
surface reflectivity (SA>0.30). While the negative ARF- 
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Fig. 9. Means of aerosol radiative forcing during the total period 
(January 2013), pollution period (daily PM 2.5 > 75pgm -3 ), and 
non-pollution period (daily PM 2.5 < 75pgm -3 ) at each site, (a) 
ARF-BOA, (b) ARF-TOA, (c) efficiency of ARF-BOA, (d) effi- 
ciency of ARF-TOA. 


TOA efficiency under the non-polluted condition could prob- 
ably be caused by faster increasing of backscattered energy 
towards the top of atmosphere in non-pollution conditions 
than under pollution ones. This faster increasing of backscat- 
tered energy towards TOA is not affected by aerosol scatter- 
ing and absorption processes. 
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Date (January 2013) 

Fig. 10. Aerosol radiative forcing at the top of the atmosphere 
(ARF-TOA) and its 5th and 95th percentile box plots during Jan- 
uary 2013 at five sites in North China Plain. 

4 Summary 

The aerosol optical properties and the radiative forcing over 
North China Plain during the intense pollution month in Jan- 
uary 2013 can be summarised as follows. 

Aside from the strong emissions based upon natural re- 
sources, meteorological conditions constitute a major con- 
tributor to the most intense pollution occurrence in January 
2013. Relative humidity and wind are two important factors 
affecting the PM 2.5 and aerosol optical properties. 

The AOD at urban, suburban, and rural stations shows con- 
sistent variations during the pollution month in 2013, which 
indicates that serious pollution not a local but a regional phe- 
nomenon in North China Plain. AOD during the pollution 
period is 2 to 4.3 times higher than in the non-pollution pe- 


riod. Alpha was larger than 0.80 at all stations and for most 
of the time, which suggests that the small aerosol particle 
size was dominant during the 2013 pollution month. The fine 
mode AOD during the pollution period at urban sites is about 
a factor of 2.5 compared with the non-pollution period, and a 
factor of ~ 5 times the values in the suburban and rural sites. 
The fine mode fraction of AOD675nm (FMF) is larger than 
80 % at all sites, indicating that the fine mode particles were 
a main contributor to the atmospheric extinction during Jan- 
uary 2013 over North China Plain. FMF during pollution pe- 
riods is obviously larger than in non-pollution periods, which 
suggests that fine mode particles have a greater extinction 
effect during pollution periods. The absorption AOD 675 nm 
(AAOD) is only ~ 0.01 at rural sites but ~ 0.03-0.07 in sub- 
urban and urban sites during pollution periods, which regis- 
ters the effect of more absorption aerosol particles at subur- 
ban and urban sites. 

The single scattering albedo (SSA) varied between 0.87 
and 0.95 during January 2013 at North China Plain. The fact 
that higher SSAs were observed during the pollution period 
than that during the non-pollution period points to the effect 
of fine mode particle growth and thereby the increase in the 
light scattering coefficient. 

The size distribution shows an obvious tri-peak pattern 
during the most intense pollution period from 10 to 16 Jan- 
uary. The fine mode effective radius during pollution periods 
is about 0.01-0.08 pm larger than that found in non-pollution 
periods, reflecting fine mode particle hygroscopic growth. 
Fine and coarse mode particle volumes during pollution pe- 
riods are ~ 0.03-0.12 pm 3 and 0.01-0.05 pm 3 larger than in 
non-pollution periods, indicative of more aerosol loading in 
the atmosphere during pollution periods. 

The ARF-BOA varies from —32 to — 144 Wm -2 , from 
—43 to — 144 W m 2 , and from —9 to —63 W m 2 for the 
whole, pollution and non-pollution periods, respectively. The 
ARF-TOA is smaller than the ARF-BOA, the former vary- 
ing from 11 to — 54 Wm -2 , 18 to —54 W m 2 and —3 
to — 21 Wm -2 for the whole, pollution and non-pollution 
periods, respectively. During the January intense period, 
the ARF-BOA exceeded — 50 Wm -2 at rural and —180 to 
— 200 Wm -2 at suburban and urban sites. The ARF-TOA 
varied with the same trends as the ARF-BOA with the ARF 
— 30 Wm -2 at rural and — 40-60 Wm -2 at urban locations. 
Positive ARF-TOA at Huimin indicates that the aerosol par- 
ticles generated a heating effect at the top of the atmo- 
sphere. This reflects the presence of a period with a large 
surface albedo resulting from a snow cover around the station 
during January 2013. The ARF-BOA efficiencies exceeded 
— 130 Wm -2 during the whole, pollution and non-pollution 
period in North China Plain. The ARF-TOA efficiencies 
in rural sites are greater than at the suburban and urban 
sites with values of —121 zb 21 Wm -2 , — 114 d= 23 Wm -2 , 
and — 134zblOWm -2 during the whole pollution and 
non-pollution period, respectively, compared to —64 to 
-87 Wm“ 2 , -55 to -77Wm“ 2 , and -67 to -lOOWm” 2 
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at suburban and urban locations. The ARF-TOA efficiency 
at Huimin during the pollution period is about 18.89 zb 
18.85 Wm -2 compared to — 32. 36 ±40.78 Wm -2 during 
non-pollution periods. 

In this article, the aerosol optical property data could be 
processed near real time (NRT), which is available as one of 
the most important input factors of haze data assimilation. 
These NRT aerosol optical property data could also be used 
in the validation of AOD and aerosol radiative forcing simu- 
lated by haze-fog modelling because their high temporal res- 
olution and continuous measurement. The above applications 
both benefits to the improvement of haze-fog forecasting sys- 
tem and therefore accurate evaluation of radiative forcing and 
radiative feedbacks on haze weather system due to haze par- 
ticles. 
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